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ABSTRACT 

During observations with the Rossi X-ray Timing Explorer from June 31st to July 
3rd 1997 we discovered two simultaneous kHz QPOs near 500 Hz and 860 Hz in the 
low-mass X-ray binary and Z source Cygnus X-2. In the X-ray color-color diagram and 
hardness-intensity diagram (HID) a clear Z track was traced out, which shifted in the 
HID within one day to higher count rates at the end of the observation. Z track shifts 
are well known to occur in Cygnus X-2; our observation for the first time catches the 
source in the act. A single kHz QPO peak was detected at the left end of the horizontal 
branch (HB) of the Z track, with a frequency of 731±20 Hz and an amplitude of 4.7to;6% 
rms in the energy band 5.0-60 keV. Further to the right on the HB, at somewhat higher 
count rates, an additional peak at 532±43 Hz was detected with an rms amplitude of 
3.0lJ;?%. When the source moved down the HB, thus when the inferred mass accretion 
rate increased, the frequency of the higher-frequency QPO increased to 839±13 Hz, and 
its amplitude decreased to 3.5 + n'o% rms. The higher-frequency QPO was also detected 
on the upper normal branch (NB) with an rms amplitude of 1.8j™% and a frequency of 
1007±15 Hz; its peak width did not show a clear correlation with inferred mass accretion 
rate. The lower-frequency QPO was most of the time undetectable, with typical upper 
limits of 2% rms, no conclusion on how this QPO behaved with mass accretion rate can 
be drawn. If the peak separation between the QPOs is the neutron star spin frequency 
(as required in some beat-frequency models) then the neutron star spin period is 2.9±0.2 
ms (346±29 Hz). This discovery makes Cygnus X-2 the fourth Z source which displays 
kHz QPOs. The properties of the kHz QPOs in Cygnus X-2 are similar to those of other 
Z sources. 

Simultaneous with the kHz QPOs, the well known horizontal branch QPOs (HBOs) 
were visible in the power spectra. At the left end of the HB the second harmonic of the 
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HBO was also detected. We also detected six small X-ray bursts. No periodic oscillations 
or QPOs were detected in any of them, with typical upper limits of 6-8 % rms. 

Subject headings: accretion, accretion disks — stars: individual (Cygnus X-2) — stars: 
neutron — X-rays: stars 
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1. Introduction 

The low-mass X-ray binary (LMXB) and Z source 
(Hasinger & van der Klis 1989) Cygnus X-2 is one 
of the best studied LMXBs. In the X-ray color-color 
diagram (CD) a Z shaped track is traced out. Its 
branches are called horizontal branch (HB), normal 
branch (NB), and flaring branch (FB). The rapid X- 
ray variability in Z sources is closely related to the 
position of the sources on the Z. Quasi-periodic oscil- 
lations (QPOs) are present on the HB and the upper 
NB with a frequency between 18 and 60 Hz (called 
HBO). QPOs with a frequency between 5 and 20 Hz 
are present on the lower part of the NB and on the be- 
ginning of the FB (called N/FBOs). Both the X-ray 
spectral changes and the rapid X-ray variability are 
thought to be due to changes in the mass accretion 
rate (M) (e.g. Hasinger & van der Klis 1989; Lamb 
1991), which is inferred to be lowest on the HB, in- 
creasing along the Z and reaching the Eddington mass 
accretion rate at the NB/FB vertex. 

Recent studies of Cygnus X-2 (EXOSAT: Kuulk- 
ers et al. 1996; Ginga: Wijnands et al. 1997a) show 
that Cygnus X-2 displays motion of the Z track in 
the hardness- intensity diagram (HID). Following Ku- 
ulkers et al. (1996), we call the different positions 
of the Z track in the HID different "overall intensity 
levels" . When the overall intensity level changes, the 
morphology of the Z changes and the whole Z moves 
in the CD as well as in the HID. Wijnands et al. 
(1997a) found that the X-ray timing behavior on the 
lower part of the NB differs significantly for different 
overall intensity levels. The fact that the X-ray spec- 
tral and X-ray timing behavior change with changing 
overall intensity level suggests that something else in 
addition to M determines the Z track and the tim- 
ing properties of Cygnus X-2. A precessing accretion 
disk, strongly suggested by the 78 day period found in 
Cygnus X-2 using the ASM/RXTE long-term X-ray 
light curve (Wijnands, Kuulkers, & Smale 1996), is a 
possible explanation of the long-term changes. 

In numerous LMXBs, mostly atoll sources (Hasinger 
& van der Klis 1989), kHz QPOs have been observed 
(see van der Klis 1997 for a recent review on kHz 
QPO properties). So far, kHz QPOs have been found 
in three Z sources (Sco X-l: van der Klis et al. 1996a, 
1997; GX5-1: van der Klis et al. 1996b; GX17+2: 
Wijnands et al. 1997b). In this letter, we report the 
discovery of two simultaneous kHz QPOs in the Z 
source Cygnus X-2. 



2. Observations and Analysis 

We observed Cygnus X-2 with the Rossi X-ray 
Timing Explorer (RXTE) from June 30th until July 
3th 1997. A total of 108 ksec of data was obtained. 
Data were collected simultaneously with 16 second 
time resolution in 129 photon energy channels (effec- 
tive energy range 2-60 keV), and with 122 (is time 
resolution in two channels (2-5.0 keV, 5.0-60 keV). 
A CD and HID were made using the 16 s data, and 
power density spectra were calculated from the 122 
,us data using 16 s data segments. 

For measuring the kHz QPO properties we fitted 
the 128-4096 Hz power spectra with a constant plus 
a broad sinusoid representing the deadtime-modified 
Poisson level (Zhang 1995, Zhang et al. 1995), and 
one or two Lorentzian peaks representing the kHz 
QPOs. For measuring the HBO we fitted the 8- 
256 Hz power spectra with a constant plus a power 
law representing the continuum, and one or two 
Lorentzian peaks representing the HBO and its sec- 
ond harmonic. Differential dead time effects (van der 
Klis 1989) were negligible due to the relatively low 
count rates. The errors on the fit parameters were 
determined using A\ 2 — 1.0; upper limits were deter- 
mined using A\ 2 = 2.71, corresponding to 95% con- 
fidence levels. The kHz QPOs were only detected in 
the 5.0-60 keV band. In the 2-5.0 keV and the com- 
bined 2-60 keV energy bands the kHz QPOs were un- 
detectable. Therefore, we used the 5.0-60 keV band 
throughout our analysis. 

For the CD we used for the soft color the logarithm 
of the count rate ratio between 3.5-6.4 keV and 2.0- 
3.5 keV, and for the hard color the logarithm of the 
ratio between 9.7-16.0 keV and 6.4-9.7 keV. For the 
HID we used as intensity the logarithm of the count 
rates in the energy band 2.0-16.0 keV, and for the 
hardness we used the hard color from the CD. For 
^12% of the time detector 4 or 5 was off. For the CD 
and HID we only used the data of the three detectors 
which were always on in order to get the position of 
Cygnus X-2 on the Z during the whole observation. 
For the power density spectra we used all available 
data. 

We used the S z parameterization (see Wijnands 
et al. 1997a and references therein) for measuring the 
position along the Z (HB/NB vertex at S z =l, NB/FB 
vertex at S z =2). Due to the fact that the position on 
the Z could be much better determined in the HID 
than in in the CD we used the S z parameterization 
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for the Z track in the HID. Thus the S z values quoted 
in this paper represent the position of the source on 
the Z track in the HID. By using logarithmic values 
for the colors and the intensity, S z does not depend 
on the values of those quantities but only on their 
variations (Wijnands et al. 1997a). 

3. Results 

In Fig. [I] the CD and the HID are shown. Cygnus 
X-2 traced out an almost complete Z track. The HID 
(Fig. |l]b) clearly shows that part of the data near 
the HB/NB vertex (regions 8, 9 and 11 in Fig. (j]b) 
was shifted to higher count rate with respect to the 
rest. No sudden jump in count rate is visible in the 
X-ray light curve so Cygnus X-2 moved smoothly to 
higher count rates. At the end of the observation the 
source moved in 0.7 days from the HB/NB vertex to 
the left of the HB and back to the vertex again which 
by then had shifted to higher count rates. This is the 
first time we catch the source in the act. Previously, 
shifted Z's were only detected on time scales of weeks 
to months (between observations), not during one ob- 
servation. In the CD (Fig. |l]a) the shift is not visible 
due to the fact that this shift takes place parallel to 
the NB (see also Wijnands et al. 1997a) and there- 
fore all data overlap. Also, the error bars in the CD 
are large, giving larger uncertainties than in the HID. 
Some of the timing properties have previously been 
observed to change when the Z track shifts in the CD 
and HID (Wijnands et al. 1997a). For the small shift 
in the Z track in our observation these changes were 
undetectable for the QPOs discussed in this paper. 

Between 2 July 1997 03:36 and 08:19 UT (5.8 ksec 
exposure time) we detected two simultaneous kHz 
QPOs (Fig.||a) with frequencies of 516±27 Hz and 
862±11 Hz (peak separation of 346±29 Hz), at 3.4ct 
and 4.4cr significance respectively. The FWHM and 
rms amplitude were 170t^ Hz and 3.6±g;|% for the 
lower-frequency QPO, and 94±H Hz and 3.5±0.4% 
for the higher- frequency QPO. Upper limits to the 
rms amplitude of the QPOs for the same time in- 
terval in the energy range 2-5.0 keV were 2.4% and 
1.1% rms, respectively. Simultaneously with these 
kHz QPO we detected the HBO. Between 2 July 1997 
09:37 and 12:57 UT (7 ksec exposure time) we de- 
tected a single kHz QPO (Fig. ||b) with a frequency 
of 779±16 Hz at 4.9er significance, with a FWHM of 
177±H Hz ' and an amplitude of 4.7±gJ% rms. The 
upper limit on the amplitude of the QPO in the 2-5.0 



keV band was 2.7% rms. kHz QPOs were detected 
at several other times, although usually below a 3<r 
significance level. 

The Z track in the HID was not homogeneously 
covered. In order to get at all positions on the Z 
track enough statistics to detect the kHz QPOs or de- 
termine useful upper limits we selected power spectra 
in the regions of the HID indicated in Fig. |l]b. Af- 
terwards we determined the average S z value of the 
power spectra selected. The frequency and the rms 
amplitudes of the kHz QPOs versus S z are shown in 
Fig. |3[ At the left most end of the HB (lowest count 
rates) only the higher-frequency QPO could be de- 
tected. When the source moved to somewhat higher 
count rates a second, lower-frequency kHz QPO be- 
came detectable. Further to the right on the HB 
the lower-frequency kHz QPO was undetectable again 
with upper limits of 2-3%. Near the HB/NB vertex 
the higher-frequency QPO was detected, but further 
down the NB the higher-frequency peak became un- 
detectable, with upper limits of 2-4% rms. With in- 
creasing M the frequencies of the higher-frequency 
kHz QPOs increased and its amplitude decreased 
(Fig. |a and c). The FWHM of the peak did not 
show a clear correlation with M . Due to the small 
number of detections of the lower-frequency QPO no 
conclusions can be drawn about the behavior of this 
QPO with M, although there are indications that its 
frequency also increases with M (see Fig. ||a). It was 
not possible to determine if the kHz QPO properties 
changed significantly as the Z track moved in the HID 
(regions 6, 7, and 12 versus 8 and 9 in Fig. [ij}). This 
was also true for the properties of the HBO. 

Simultaneously with the kHz QPOs we detected 
the HBO. Its properties (in the 5.0-60 keV band) 
versus S z are also shown in Fig. ||. The HBO 
second harmonic was only detectable in region 1 
(S z =-0.09±0.1) of Fig. [l]b with a frequency of 73.6±0.8 
Hz, an rms amplitude of 3.1+^%, and a FWHM 
of 13±3 Hz. The HBO fundamental in this region 
had a frequency of 36.2±0.2 Hz, an rms amplitude 
of 4.5t^%, and a FWHM of 6.8±%$ Hz. When the 
source moved further to the right on the HB the sec- 
ond harmonic became undetectable with typical up- 
per limits on the amplitude of 1-2% rms. The HBO 
fundamental was detected down to about halfway 
the NB (region 14; S z =1.46±0.05). Its frequency in- 
creases from 36.2±0.2 Hz at the left end of the HB to 
56.4±0.1 Hz at the HB/NB vertex (Fig. |b). After 
that its frequency remained approximately constant. 
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The rms amplitude of the HBO first decreased on the 
HB with increasing M, but around S z =0.6 it started 
to increase again (Fig. ||d). At around S z =l.l it de- 
creased again with increasing M until it became un- 
detectable further down the NB. Its FWHM also first 
decreased on the HB with increasing M (Fig. ||f), and 
at the same position on the Z where the amplitude of 
the HBO started to increase again with increasing M, 
the FWHM also started to increase again. 

We detected six bursts in the X-ray light curve 
(typical durations of 4-6 s) , which are very similar to 
those found in EXOSAT (Kuulkers, van der Klis, & 
van Paradijs 1995) and Ginga (Wijnands et al. 1997a) 
data. We searched for periodic oscillations and QPOs 
during the bursts, but none were detected (typical 
upper limits of 6-8% rms). A more detailed study of 
the X-ray bursts will be reported elsewhere. 

4. Discussion 

We have discovered two simultaneous kHz QPOs 
near 500 Hz and 860 Hz in the Z source Cygnus X-2. 
The frequency of the higher-frequency QPO increased 
when the source moved from the left end of the HB to 
the upper NB, and there are indications that the same 
is true for the frequency of the lower- frequency QPO. 
Cygnus X-2 is the fourth Z source which displays kHz 
QPOs. Most likely the remaining two Z sources (GX 
340+0 and GX 349+2) will also display kHz QPOs 
when they are observed on the (left end of the) HB 
(note that GX 349+2 has never been observed in the 
HB). 

The properties of the kHz QPOs in Z sources are 
very similar to each other. In each source the kHz 
QPOs increase in frequency and the higher-frequency 
QPO decreases in amplitude when the sources move 
down the Z, thus when M increases. The amplitude 
of the lower-frequency QPO and the FWHM of the 
peaks do not have a clear correlation with M. The 
kHz QPOs in Z sources are also very similar to what 
is observed (e.g. Strohmayer et al. 1996; Zhang et al. 
1996; Berger et al. 1996; Smale et al. 1997) in other, 
less luminous LMXBs (the atoll sources [Hasinger & 
van der Klis 1989]). In both type of sources the fre- 
quency of the kHz QPOs increases with inferred M, 
the kHz QPOs are strongest in the highest photon en- 
ergy band, the maximum frequencies so far detected 
are between 1000 and 1200 Hz, and the separation 
between the kHz QPOs lies between 250 and 370 Hz 
(see van der Klis 1997 for a recent review on kHz QPO 



properties). Therefore, it seems likely that kHz QPOs 
in atoll sources and in Z sources are produced by the 
same mechanism. 

In Sco X-l (van der Klis et al. 1996a; 1997) the 
peak separation between the two kHz QPOs decreases 
with increasing M[]. In GX 17+2 a similar decrease 
in peak separation could not be excluded (Wijnands 
et al. 1997b) and due to the sparcety of detections 
with two simultaneous kHz QPO in the present data 
nothing can be said for Cygnus X-2. If the peak sep- 
aration was a measure for the neutron spin frequency 
then the spin frequency in Cygnus X-2 is 346±29 Hz 
(spin period of 2.9+0.2 ms). In the four Z sources 
for which kHz QPOs have been discovered the kHz 
QPOs were detected simultaneously with the HBO, 
ruling out the possibility that both the HBOs and 
the kHz QPOs are magnetospheric beat frequencies 
(van der Klis et al. 1997b). 

At the left end of the HB the amplitude, in the en- 
ergy range 5-60 keV, of the higher-frequency QPO in 
Cyg X-2 is 4.7+o;g% rms. This amplitude is smaller 
than the typical amplitudes of the kHz QPOs in the 
atoll sources, which are thought to have significantly 
weaker magnetic fields, as predicted by the sonic- 
point model (Miller, Lamb, & Psaltis 1997). The 
QPOs in Cyg X-2 are similar in strength to those in 
GX 17+2 (2-5% rms in 5-60 keV), even though GX 
17+2 may have a weaker magnetic field, as suggested 
by its weaker HBO and by X-ray spectral modeling 
(Psaltis, Lamb, & Miller 1995). This indicates that 
if the strength of the stellar magnetic field affects the 
amplitudes of the kHz QPOs, it is not the only vari- 
able that does so, and that other variables, such as the 
multipolar structure and orientation of the magnetic 
field, and the spin rate of the star, also play a role. If 
the orientation of the magnetic field with respect to 
the plane of the disk plays a role its affect should be 
visible in the amplitude of the kHz QPOs in Cygnus 
X-2 when these QPOs are detected during different 
overall intensity levels, thus at different phases of the 
78 day long-term X-ray cycle. 

Note added in manuscript After submission of this 
paper we discovered also two simultaneous kHz QPOs 
in the Z source GX 340+0 at 348±^ Hz and 722+13 
Hz (peak separation of 374+24 Hz, when the source 
was on the horizontal branch (Jonker et al. 1997). 
The FWHM and rms amplitude (in the energy range 

Recently, Mendez et al. (1997) found that also in the atoll 
source 4U 1608—52 the peak separation may not be constant. 
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5.0-60 keV) were 114^ Hz and 2.5±g;^% for the 
lower-frequency QPO, and 177±£| Hz and 3.8±g;t% 
for the higher-frequency QPO. 
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0.36 0.38 0.4 0.42 
Soft Color 

Fig. 1. — Color-color diagram (a) and hardness- 
intensity diagram (b) of Cygnus X-2. The soft color 
is the logarithm of the count rate ratio between 3.5- 
6.4 keV and 2.0-3.5 keV; the hard color is the log- 
arithm of the count rate ratio between 9.7-16.0 keV 
and 6.4-9.7 keV; the intensity is the logarithm of the 
3-detector count rate in the photon energy range 2.0- 
16.0 keV. Both diagrams are corrected for background 
(~50 cts/s in the energy range 2.0-16.0 keV); the 
count rate is not dead-time corrected (4-6% correc- 
tion). All points are 64 s averages. Typical error 
bars are shown at the bottom right corners of the di- 
agrams. The boxes in the HID are the regions which 
were selected to study the timing behavior. 
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Fig. 2. — Typical Leahy normalized power spectra in 
the energy range 5.0-60 keV showing the kHz QPOs. 
The upwards slope at kHz frequencies is due to instru- 
mental dead time, (a) is the power spectrum for the 
interval 2 July 1997 03:36-08:19 UT; (b) the power 
spectrum for the interval 2 July 1997 09:37-12:57 UT. 



9 




Fig. 3.— (a) Frequency of the kHz QPOs, (b) fre- 
quency of the HBO fundamental, (c) the rms ampli- 
tude of the higher- frequency kHz QPO, (d) the rms 
amplitude of the HBO fundamental, (e) the rms am- 
plitude of the lower-frequency kHz QPO, and (f) the 
FWHM of the HBO fundamental as a function of 
S z . In (a) the squares represent the higher- frequency 
QPO and the circles the lower-frequency QPOs. Open 
symbols represent detections of the kHz QPOs with 
significance between a 2 and 3a, the filled points rep- 
resent detections with >3cr significance. 
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